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In two recent articles |11,[2|], Meiser and Meystre describe the coupled dynamics of a dense gas of atoms and 
an optical cavity pumped by a laser field. They make two important simplifying assumptions: (i) the gas of 
atoms forms a regular lattice and can be replaced by a fictitious mirror, and (ii) the atoms .strive to minimize the 
dipole potential. We show that the two assumptions are inconsistent: the configuration of atoms minimizing the 
dipole potential is not a perfect lattice. Assumption (ii) is erroneous, as in the strong coupling regime the dipole 
force does not arise from the dipole potential. The real steady state, where the dipole forces vanish, is indeed a 
regular lattice. Furthermore, the bistability predicted in |T, 2] does not occur in this system. 
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In two recent articles f?, '5], Meiser and Meystre describe 
the coupled dynamics of movable atoms, mirrors and light in a 
pumped optical cavity. They take the effect of the atoms on the 
light field into account using a 1 -dimensional non-perturbative 
model introduced by Deutsch et al. |3]. The mechanical ef- 
fects of the off-resonant light on the atoms are described by 
a dipole Hamiltonian. They claim that the atoms trapped in 
the cavity field self-organize to form a regular lattice of atom 
clouds which behaves as a beam splitter (BS). This BS ef- 
fectively splits the single cavity into two coupled resonators 
("left" and "right"). An important prediction of iH 0] is a 
bistability effect: optical forces will push this "atom BS" to a 
position where it is approximately an integer multiple of the 
half-wavelength away from the left (right) cavity mirror. Thus 
for a certain parameter regime the "left" ("right") cavity is on 
resonance with the pump laser, and has intense light, while the 
other cavity has weak field. 

We argue that the model used by Meiser and Meystre con- 
tains contradictions and errors involving the way the dynam- 
ics of the atoms (taking place on a much shorter timescale 
than that of the movable cavity mirror in |[l|]) is treated. The 
model hinges on two key assumptions. The first is a claim 
based on |3]: (i) the gas of atoms forms a regular lattice of 
pancake-shaped clouds with lattice constant do = A/2(l + 
2atanA/7r), and can thus be replaced by a fictitious "atom 
BS". Here A — kria/{2€o) is the dimensionless polarizabil- 
ity density of a cloud of surface density 77 composed of atoms 
of polarizability a < 0, with k = 2tt/X denoting the free- 
space wavenumber of the pump laser The second assumption 
is used to find the steady state position of the atom BS: (ii) 
the atoms strive to minimize the total dipole potential Hmt 
(eq. (2) of [3], or eq. (3) of [2]). In this Comment we show 
that these assumptions are inconsistent, and that (ii) has to be 
replaced by a formula for the force on the atoms. We present 
the required expression for the force which allows us to deter- 
mine the steady state of this system. This does not exhibit the 
bistability phenomenon of 1 1, 2]. 

The setup considered in 1 1 , 2] is an open system: laser light 
enters the cavity through the mirrors, carrying momentum and 
energy, and is coupled out of the cavity at the mirrors. It 
is not at all clear what function of the system parameters is 
minimized in a steady state. Minimizing the dipole energy 
certainly does not lead to steady state configurations. To il- 
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FIG. 1: (color online) A configuration found by Monte-Carlo type 
minimization of the dipole potential, for A'^ = 100 atom clouds with 
polarizability A = 0.1, symmetrically pumped in free space (no cav- 
ity). The local lattice constant (blue dots, in units of the original 
do ~ 0.468A) and the electric field intensity at the atomic posi- 
tions (red line, exponential scale, in units of the pump intensity Jo) 
are plotted. The total dipole energy is some 1000 times lower than 
that of the regular optical lattice of assumption (i). The atoms form 
two slabs with lattice constant exceeding do, impervious to the pump 
laser (pump frequency deep in the band gap, see |3]). These two self- 
organized atom mirrors constitute a high-Q cavity, trapping light in 
the middle. There the intensity is some 10* times that of the pump 
lasers. Inside the "atom mirrors" the intensity falls off exponentially. 



lustrate this point, it is worthwhile to consider a conceptually 
simpler situation, namely atoms trapped in a standing wave 
laser field without any cavities involved. As illustrated in 
Fig. [T] using the model of Meiser and Meystre but allowing 
all of the atom clouds to move independently, a Monte Carlo 
algorithm to minimize the dipole energy leads to configura- 
tions differing from a simple lattice. Starting from the regular 
lattice configuration, the atoms decrease the dipole energy by 
forming a self-organized cavity resonant with the pump beam. 
The intensity inside this cavity is thus enhanced by the reso- 
nance, and a few atoms coupled to this intense field contribute 
to the total dipole potential by such a large amount, that its 
absolute value exceeds that of the original energy by orders 
of magnitude (see the figure caption for details). The mir- 
rors of this self-organized cavity are slabs of an atomic lattice 
with lattice constant exceeding do such that the pump field is 
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in the band gap fs*] of these slabs and decays exponentially 
inside them. The complicated spatial structure of this solu- 
tion clearly shows that assumption (ii) cannot lead to (i), one 
of them has to be dropped. The configuration represented in 
Fig.[T]is not a steady state at all: as we discuss below, radia- 
tion pressure would push the two slabs apart. Thus, assump- 
tion (ii) has to be revisited. We remark that inside an optical 
cavity, dipole energy minimization leads to similar artifacts, 
in the "superstrong coupling" limit as well. 

The reason why dipole energy minimization does not sup- 
ply the steady states is that the atomic positions are coupled 
parametrically to the Ught field. This problem is met, and is 
tackled in a very neat way, when the motion of atoms coupled 
to a single-mode high-Q optical resonator is to be described in 
a classical approximation (Ehrenfest theorem). In the stan- 
dard approach (e.g., |4J), the starting point is the quantum 
Hamiltonian H = X^^i i?atom , CTj , ct] ) + i?ficid(a, a^) + 

hY^^^^ 9{zj){(^]a + o'j-a^), with Zj, pj, and aj denoting the 
position, momentum, and deexcitation operator of the j-th 
atom, g{z) the mode function of the cavity, and d the cav- 
ity photon annihilation operator. The coupling between atoms 
and the cavity field is of the celebrated Jaynes-Cummings 
type. The classical approximation should furnish equations of 
motion for — {zj) and pj — {pj). In order to derive these, 
for slowly moving atoms, the separation of the timescales 
is invoked. The internal variables d and aj, which equili- 
brate fast on the timescale of atomic motion, are replaced 
by their adiabatic steady state expectation values, which de- 
pend on all of the atom coordinates z; (as well as the in- 
tensities and phases of the pumping lasers). One could be 
tempted to use the "effective" Hamiltonian Hcs{zj,Pj) ob- 
tained from H in this way, and derive the atomic dynamics 
from it via the Heisenberg equations of motion, i.e., dpj/dt = 
{—dHcff /dzj) = —d{Hos) / dzj = —d{H)/dzj. The correct 
procedure, however, is to apply the adiabatic approximation 
to the original Heisenberg equations 
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where the force operator is Fj = —h{a^jd+d-jd'^)dg{zj)/dzj, 

and use dpj/dt = {Fj). In other words, the differentia- 
tion should only be applied w.r.t. explicit Zj -dependence of 
(H), since (dH/dzj) is not the same as d{H)/dzj. Using 
d{H) I dzj to define the dynamics supplies steady states where 
the "dipole potential" (H) — (Hcff) is minimized, but these 
are not the true steady states of the system: in these states the 
force on the atom {Fj ) does not vanish. 

We now turn to the setup considered by Meiser and 
Meystre, where the back-action of the atoms on light is so 
substantial that the cavity no longer has a fixed mode function. 
To be self-contained, and to fix notation, we briefly summa- 
rize the model, detailed in llj, 0, [H . 

Assuming that the atoms are fixed on the timescale of the 
field dynamics, the light field is calculated by solving the 



Helmholtz equation in a one-dimensional approximation, 

a^£:(z) + k^E{z) ^ -2kKE{z) ^ 5{z - Zj). (2) 

i 

The right-hand-side embodies the polarizability of the trapped 
atoms, which are assumed to form pancake-shaped clouds of 
axial size much smaller than a wavelength. The solution of 
this equation is trivial: between two atom clouds, the electric 
field is a superposition 

E{zj-i <z< Zj) = Aje'''^'^'^^ + Bje-'^'^'-''''^ 

= Q.ie^'^^^-^^-i) + i:)j_ie"''=(^"^^-i). (3) 

The field has to fulfil boundary conditions: 

E{z = Zj - 0) = E{z = Zj + 0); (4a) 
d^E{z = Zj - 0) = d^E{z ^ z.j + 0) + 2kKE{zj). (4b) 

These conditions are equivalent to representing the atom 



clouds by BS's, i.e., Aj = tCj + rBj, Dj = tBj 



with complex reflection and transmission coefficients r = 
iA/{l-iA),t^ 1/(1 -iA) |3]. 

The dynamics of the atoms is given by the dipole force act- 
ing on them. Instead of minimizing a dipole potential, the true 
steady state of the system is then specified by the positions of 
all the atom clouds Zj,j = 1, . . . ,N such that the optical field 
of the cavity - the solution of (|2]i - exerts no net force on any 
of the clouds. In the following we show two ways to calculate 
this force acting on an infinitely thin atom cloud (BS). 

The force on an atom cloud is can be obtained by inte- 
grating the the force on a single atom over the whole cloud. 
A microscopic model of light-matter interaction leads to two 
types of force: the dispersive dipole force and the dissipa- 
tive scattering force p"]. This latter is often referred to as 
"radiation pressure", but following Meiser and Meystre we 
use this term to denote the mechanical effects of light in gen- 
eral. In iQl, 0] flie atom-pump detuning is assumed to be so 
large that the scattering force can be neglected, tantamount 
to assuming A G M. For linearly polarizable particles the 
dipole force time-averaged over an optical period is given fsj] 
F — -jaV |F(x)| . Calculating this force for an infinitely 
thin disk-shaped atom cloud poses a problem, as the elec- 
tric field E{z) is not differentiable at the atomic positions 
Zj. One must calculate the force on a disk of finite extent 
Zj — w . . . Zj + w, and only then take the limit w ^ 0. Since 
the electric field is polarized in the plane of the disk, there 
is no surface contribution [6], and the integral in the limit of 
vanishing width gives 
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(9,|F|'(zj-0) + 9,|F|'(zj+0)) (5) 



for the force on a unit surface ("radiation pressure"). This 
result is independent of the way in which the limit is ap- 
proached, i.e. of the axial density distribution of the cloud. 
Substituting the modal decomposition of Eq. Q, using the BS 
relations, and the fact that as A £ M, we have \Aj\ +\Dj\ = 



\B, 



\Cj\ , some algebra leads to the simple formula 
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There is another, macroscopic way to arrive at the light- 
induced force on a scatterer. This consists of calculating the 
Maxwell stress tensor and integrating it on an arbitrary ficti- 
tious surface enclosing the body (see, e.g., Ref. |7]). In the ID 
model of 1 1, 3], this is very easily done. For a selected atom 
cloud, we take the surface around it to consist of two planes 
orthogonal to the cavity axis, between the atom cloud and the 
two neighbouring clouds. As the electromagnetic wave in- 
side the cavity is transverse, both the E and B vectors lie in 
the planes, and the only part of the stress tensor contributing 
to the integral is the term with the energy density. For the 
plane waves of (O this results in (|6]l. We remark that this line 
of thought is also alluded to by Meiser and Meystre, and al- 
though is not applied to the atoms, it is used to derive the force 
on the movable cavity mirror in Eq. (10) of |T{]. 

The difference between (a) minimizing the dipole poten- 
tial, and (b) requiring the dipole force to vanish is illustrated 
in Fig. |2] Here we put a weakly reflective BS ( A = 0.1) in 
a high-Q cavity composed of two Dirac-(5 distributions of po- 
larizability with A = 10, corresponding to transmission prob- 
ability T w 0.01, as in f?, '2]. To a good approximation the 
cavity only supports a single sin(z)-mode, and since the BS 
is weakly reflective, we are not in the "superstrong coupling" 
limit: the Hamiltonian approach would constitute a reason- 
able approximation. Note however, that in Fig. |2ja), the in- 
tensities to the right of the trapped atom are slightly higher 
than to its left, indicating the corrections to the single-mode 
approximation, and also the bistability of 1 1]. When minimiz- 
ing the dipole energy (a), the BS occupies a position where 
it is not coupled strongly to the cavity, far from an antin- 
ode of the mode function. Thus the system is near-resonant 
with the driving field, and the intracavity intensity is enhanced 
by the resonance, some 500 times the intensity outside. Ex- 
tremely high intensity means large negative dipole energy for 
the trapped atom - however, note that at the position of the 
atom, the derivative of the intensity is nonvanishing, and thus 
the force F is nonzero. Requiring this force to vanish brings 
the atom to a mode function antinode (b), where the coupling 
is stronger, hence the frequency shift is larger, than in (a). For 
this specific example, this is already enough to shift the sys- 
tem out of resonance and to decrease the intracavity intensity 
below the free-space value. 

The BS used to model the atom cloud in the cavity consid- 
ered by fH, 121] , has an effective Awl, and thus the system is in 
the "superstrong coupling limit". We plot the force exerted by 
the cavity field on such a BS as a function of its position and of 
the detuning (equivalently, cavity size) in Fig. [3] The bound- 
aries between the gray shaded and white areas correspond to 
equilibrium. Note however, that for a fixed drive detuning, of 
the two equlibrium solutions per half wavelength only one is 
stable, the one where d/dz^ F < 0. Thus there is no bistabil- 
ity of the kind predicted in 1 1 ] . The areas enclosed by the solid 
contour lines, where the force on the BS becomes very large, 
indicate that the system is on resonance and the field gradi- 
ent at the BS position is high. These correspond to the black 
structures on Fig. 2. of [2], where the determinant D (Eq. (12) 
of ytl) becomes small, which there are falsely interpreted as 
equilibrium positions. 
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FIG. 2: (color online) Steady states of a single atom (atom cloud) 
of polarizability A — 0.1 in a symmetrically pumped cavity, ac- 
cording to (a) minimization of the dipole energy; and (b) vanishing 
dipole force. The cavity is formed by two highly reflective mirrors 
(A — 10) at 2; = and z — 1.501 A. The position of the atom 
(black dot and vertical dotted line) induces no substantial change of 
the mode function of the cavity. However, it influences the inten- 
sity of the cavity mode (filled red curve, in units of the free-space 
intensity Jo, mirrored for better visibility). In (a), the system is on 
resonance, the intracavity intensity is so high (some 500/o) that the 
intensity outside the cavity is hardly seen on this range; in (b), the 
atom is maximally coupled to the cavity, shifts it out of resonance 
and reduces the intensity below the free-space value. 



In the example shown in Fig.|2l the maxima of the intensity 
to the left and to the right of the trapped BS are equal. In fact, 
an analogous statement holds for the steady state of any one- 
dimensional system composed of consecutive beam split- 
ters, held together by the dipole force, regardless of the BS 
parameters of the system components. This follows from for- 
mula (|6]l for the dipole force, whereby for every j — 1, . . . ,N: 



.f' + \Bf' 
no absorption, \Aj\ + \Dj 
These relations imply that 

|Ci| = IC2I \Cn\; 
\Di\ = \D2\ =,...,= \Dn 



Dj\ . Furthermore, since there is 



\Bjr + \Cj\ 



at every BS. 
^il = IA2I =,...,= \An\ - 
; jSil = IS2I -,...,= \Bn\ = 
i.e., the plane waves pass the 



atom clouds unattenuated, suffering only phase shifts. Thus, 
the envelope of the intensity oscillations of the electric field 
is constant throughout the sample. This is in line with the in- 
tuitive picture of "radiation pressure" caused by the collisions 
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FIG. 3: Force on a single beam splitter with A = 1 in a standing- 
wave cavity with parameters as in Qj, Q] (mirror transmission prob- 
ability r 0) pumped by a laser via one of the end mirrors. Con- 
tour lines are: solid, |F| = lOFo; dashed, \F\ = Fo/10; dotted, 
|F| = Fo/1000, where Fo is the radiation pressure force that would 
act on the beam splitter in the absence of the cavity. Gray back- 
ground indicates positive (rightward) forces, white background nega- 
tive (leftward) forces. Extremely high forces (|F| > lOFo, the small 
areas enclosed by the solid contour lines) occur whenever a part of 
the cavity is on resonance with the drive. The boundaries between 
the gray and white areas, where f ~ 0, are the (stable or unstable) 
equilibrium positions. 



of photons with the atom clouds. However, these results are in 
direct contradiction to those obtained by Meiser and Meystre: 
the nonconstant field envelopes in the "bistability regime" are 
explicitly plotted in Fig. 3 of |2]. 

We now revisit assumption (i). Assuming that there exists 
a steady state of N identical disk-shaped atom clouds trapped 
by the light field in a cavity, we find that these clouds have to 



form a perfect lattice. This is true because in the steady state 
the light permeates the stack of clouds unattenuated. Thus 
\E{x)f = \Eof + \Eif + 2 \EoEi\cos{2kx~^{x)) every- 
where in the sample, the clouds only contribute to the phase: 
^{xj < X < Xj+i) = X/- The phase slip at the 

rth cloud xi depends not only on the polarizability density 
A of the clouds, but also on the pump asymmetry, i.e., the 
ratio of the intensities of the left- and rightwards propagat- 
ing waves, see 1^ for details. Since the clouds are identical, 
and the light fills the structure unattenuated, both these pa- 
rameters are equal for all clouds. Thus, xi = X for every 
I ^ 1, . . . ,N, therefore the atom clouds form a perfect lattice 
(possibly with gaps of an integer multiple of A/2) with lattice 
constant d = ^(^r — x)- However, the value of the phase slip 
X, and thus of the lattice constant d, is not trivial to determine, 
due to the dependence on the pump asymmetry. For the atoms 
trapped inside the asymmetrically pumped cavity considered 
by Meiser and Meystre, this asymmetry varies with the pump 
detuning: on resonance, it is negligible, while far from reso- 
nance, it is substantial. Thus Eqs. (24) and (25) of Ol] (taken 
from fs*] for symmetric pumping) cannot be applied to this 
system. Even more crucially, whether or not the equilibrium 
sets in depends on the dynamics. We have found [St] that for 
large lattices, even a small pump asymmetry can lead to a dy- 
namical instability of the equilibrium configuration. 

Common wisdom holds that "dipole force is conserva- 
tive". Even in the strong coupling regime of cavity QED 
it is often possible to construct a "potential" by integrat- 
ing the dipole force (e.g., |4]). For many particles trapped 
in the same cavity, this potential should also include the 
field-mediated (parametric) interaction between those parti- 
cles. We have found J^, that in the generic case even 
this approach breaks down. For asymmetric pumping (as in 
li i), 5Fj(zi, . . . , ZM)/dzi ^ dFiizi, ZN)/dz„ and 
thus no potential function V{zi, . . . , zm) can be constructed 
that obeys the Young theorem about the commutativity of par- 
tial derivatives. 
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